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•7^. Abstract 

^5 . The ratios of particle densities in lead-lead and proton-lead collisions to 

particle density in proton-proton collision in the central rapidity region at the 
LHC energy are predicted on the basis of wounded quark-diquark model. 
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The wounded quark-diquark model [1, 2] proved to be rather successful in description 
of particles production from nuclear targets. Assuming that high energy interac- 
tions of nucleons are dominated by independent interactions of its two constituents, 
a quark and a diquark, it was possible to describe pp, dAu, CuCu and AuAu mul- 
tiplicity data collected at the RHIC collider [3]. It indicates that in all hadronic 
collisions the early stage of the particle production process can be understood as a 
simple superposition of contributions from hadronic constituents. As explained in 
[2] this does not preclude further collective evolution of the system that is obviously 
present [4, 5]. 

Encouraged by these results we present here quantitative predictions of the 
wounded quark-diquark model for the particle density ratios Rab = NAs/Npp in the 
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central rapidity region of PbPb and pPb collisions at the LHC energy y^ = 5500 
GeV.i 

Our main conclusion is that the model provides rather precise predictions for 
the nuclear collisions at LHC energies. This should allow its effective test when the 
data are available. 

In the next section the prediction of the wounded quark-diquark model for par- 
ticle density in the central rapidity region in PbPb collision is presented. In section 
3 we discuss the consequences of the model for midrapidity density in pPb colli- 
sions. Our conclusions are listed in the last section where also some comments are 
included. 

2 PbPb collision 

The relation between particle production in nucleon-nucleon and symmetric nucleus- 
nucleus collisions implied by the wounded quark-diquark model is given by [1] 

P _ NphPbivl _ n/g^ 

where the R.H.S of this equation is independent of rapidity y.^ Npbpt,{y) and Npp{y) 
are the particle densities in PbPb and pp collisions, respectively. Wp is the average 
number of wounded constituents in a single pp collision (per one proton). Mean 
number of wounded quarks and diquarks in both colliding nuclei Wp^pl at a given 
impact parameter b is given by (mass number A = 208) [6] 

^ptplib) = ^^ / Tib -s){2-[l- pMs^ - [1 - pMs)]^} dh, (2) 

with G{s) defined as 

G{s) = f (fs'ainis - s')T{s'), (3) 

where T{s) is the nuclear thickness function T{s) = f dzp{\^s'^ + z"^) (normalized 
to unity). Here and in the following for the nuclear density p we take the standard 
Woods-Saxon formula with the nuclear radius Rpb = 6.5 fm and the skin depth 
d = 0.54 fm. apbPb{b) is the inelastic differential PbPb cross section.^ Finally, Pq 



"'^The density of particles produced in pp collision at the LHC energy cannot be predicted in the 
present approach. 

^Provided we are far enough from the fragmentation regions, where contributions from cascade 
and unwounded constituents are expected [2] . 

^o'PbPbib) = 1, except at very large impact parameters (b > 14 fm) which are of no interest. 



and pd are the probabilities for a quark and a diquark to interact in a single pp 
collision, respectively. 

We assume the differential inelastic pp cross section (Tin(s) (probability for in- 
elastic pp collision at a given impact parameter s) to be in a simple Gaussian form^ 

^in(s) = e-^'/-^ (4) 

where x^ = o"in/vr and o-jn is the total inelastic pp cross section o-jn = Jo'iii(s)(i^s. 

The multiplicity data are usually presented versus the number of wounded nu- 
cleons [8] 

4t.W = ^^ Inh - s) {l - [1 - G{s)f] dh. (5) 

o'pbPbyo) J I- ■> 

This completes all necessary formulas. 

To obtain Wp we followed exactly the procedure proposed at lower energies, 
where we extracted this number [1] by studying differential elastic pp scattering 
cross section data. Indeed, assuming a nucleon to be composed of a quark and a 
diquark, it was possible to describe the small momentum transfer, |t| < 3 GeV^, 
elastic pp and np scattering cross section data with a very high precision [9]. In 
the present case we studied the small t elastic pp scattering data at the Tevatron 
energy giving Wp = 1.24±0.01. Considering many different predictions regarding 
elastic pp scattering at 14000 GeV [10] we obtained Wp = 1.28 ± 0.02. Thus in 
our calculations at ^/s = 5500 GeV for the average number of wounded quarks and 
diquarks in a single pp collision, per one colliding proton, we take 

^«+^) = 1.26 ±0.02. (6) 

This number is the dominant uncertainty of our approach. The detailed discus- 
sion of this problem, however, is beyond the scope of this investigation. 

Since the total inelastic pp cross section o"in is not known at ^/s = 5500 GeV we 
performed our calculations for three different inelastic cross sections (Tin = 60, 67 and 
75 mb. We noticed that at a given number of wounded nucleons we obtain practically 
the same number of wounded constituents for each value of (Tin. This observation 
allows for predictions at the LHC energy, which are practically independent of the 
value of (Tin. 

The calculated numbers for ain = 60 mb are presented in Table 2, where following 
[9], we assumed Pq = Pdf^ = Wp /3.^ 

Dividing Wpf^pl by 2wp we obtain our prediction for the ratio Raa (1) shown 
in Fig. 1. For comparison the prediction of the wounded nucleon model [8] is also 
shown. 



^We believe that crin(O) = 1 is very close to reality. At ISR energies crin(O) = 0.92 [7]. 

^We also checked different choices, ranging from pd = Pq to pd = 2pq. We observe that the 
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relation Wpf^pl vs Wp^pf^ is not changed 
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Table 1 : Mean number of wounded quarks and di quarks Wp^^pl and wounded nucle- 
ons Wplp^ in PbPb collision as a function of the impact parameter b. 

3 pPb collision 



At the vanishing cm. rapidity we have the following relation between particle 
production in pp and pA collisions [1] 
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where the average number of wounded quarks and diquarks in pPb collision w p^^ 
at a fixed impact parameter b is given by 
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The first term gives the number of wounded constituents in the target {Pb nu- 
cleus). Indeed, it is the number of wounded nucleons in the target times the number 
of wounded constituents in a single pp collision. The second term gives the number 
of wounded constituents in the projectile that underwent many inelastic collisions. 
The derivation of this term is presented in the appendix. 

Mean number of wounded nucleons at a given impact parameter b is given by 
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where the first term gives the number of wounded nucleons in the target, plus one 
wounded nucleon being the projectile itself. 
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Figure 1: Wounded quark-diquark model prediction for the multiplicity ratio of 
particles produced in PbPb collision to those produced in pp collision at any rapidity 
y. The grey band reflects the uncertainty in the value of Wp . The prediction of 
the wounded nucleon model is also shown. 

Again, we performed the calculations for three different inelastic cross sections 
(Tin = 60, 67 and 75 mb. At a given impact parameter b we obtain significantly 
different numbers of wounded nucleons and wounded constituents, however, when 
we plot Wpp^ vs Wppj the three curves almost exactly follow each other. Similarly 
to the previous case of PbPb collision, this observation allows for predictions at the 
LHC energy which are independent of the value of o"in. 

The obtained numbers for ai^ = 75 mb^ and Pg = Pd/2 = Wp /3 are presented 
in Table 2. 

Dividing w^p^ by 2wp we obtain our prediction for the ratio RpA (7) pre- 
sented in Fig. 2. The maximal number of wounded nucleons is 14 and 17 for a-^ = 60 
and 75 mb, respectively. For comparison we also show the prediction of the wounded 
nucleon model. 

It is not surprising that the wounded quark-diquark model prediction is rather 
close to the line predicted by the wounded nucleon model. Indeed, comparing both 
scenarios the only difference is the projectile that undergoes many inelastic collisions 
producing slightly more particles [2]. 
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This time we take the largest number. Maximal mimbcr of wounded nucleons noticeably 
depends on (Tin. The relation w^'pf^ vs w^p;, hardly depends on it, however. 
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Table 2: Mean number of wounded quarks and diquarks w^pf^ and wounded nucle- 
ons Wppft in pPb collision as a function of the impact parameter b. 

4 Conclusions and comments 



Our conclusions can be formulated as follows. 

(i) Encouraged by a very good agreement of the wounded quark-diquark model 
with the RHIC pp, dAu, CuCu and AuAu data, we evaluated particle densities in 
the central rapidity region in PbPb and pPb collisions at the LHC energy y/s = 5500 
GeV. 

(ii) In our approach the particle density in PbPb (at the central rapidity region) 
and pPb (at midrapidity) is proportional to the density of particles produced in 
an elementary pp collision. Since the pp particle density is presently unknown and 
it cannot be calculated in the present approach we only give the ratios Rab = 

Nab/Npp. 

(iii) The dominant uncertainty of our calculation is the number of wounded 
quarks and diquarks in a single pp collision at ^/s = 5500 GeV which we estimated 
to be 1.26 ±0.02. 

(iv) Since the total inelastic pp cross section is not known at ^/s = 5500 GeV we 
performed our calculations for three different inelastic cross sections (Xin = 60, 67 and 
75 mb. The functional relation between number of wounded quarks and diquarks 
and number of wounded nucleons practically does not depend on the value of ai^. 
This observation allowed for predictions at the LHC energy, which are independent 
of the value of (Tin. 

Following comments are in order. 

(a) Our prediction regarded the multiplicity density ratio Raa can be also applied 
to the total multiplicities measured for central PbPb collisions. For such centralities 
additional contributions from cascade and unwounded constituents seem to be less 
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Figure 2: Wounded quark-diquark model prediction for the multiplicity ratio at 
midrapidity of particles produced in pPb collision to those produced in pp collision. 
The maximal number of wounded nucleons is 14 and 17 for o"in = 60 and 75 mb, 
respectively. The grey band reflects the uncertainty in the value of Wp . The 
prediction of the wounded nucleon model is also shown. 

important [2]. 

(b) We found previously that the 200 GeV RHIC data in the range \y\ < 3.7 
can be solely described by the contribution from the wounded constituents. Beyond 
this region unwounded constituents and cascade seems to appear [2]. Assuming that 
these additional contributions begin at y proportional to the rapidity beam Y, it 
allows us to estimate that at ^/s = 5500 GeV the ratio Raa should be independent 
of y in the approximate range \y\ < 6. 

(c) In principle our predictions could be applied to any energy provided Cin 
remains in the range from 60 mb to 75 mb. The only difference is the number of 
wounded quarks and diquarks in a single pp collision. For instance at ^/s = 14000 
GeV we estimate wl^^'^^ to be 1.28 ± 0.02. 
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A Appendix: Wounded constituents in the pro- 
jectile 

The average number of wounded quarks and diquarks in a nucleon that underwent 
exactly k inelastic collisions is given by 

Wk = l-{1-Pg)' + 1-{1-Pd)\ (10) 

where Pg and p^ are the probabilities for a quark and a diquark to interact in a single 
pp collision, respectively. 

The probability that the nucleon at a given impact parameter b underwent ex- 
actly k inelastic collisions is given by a standard formula 

where we assume that at least one inelastic collision takes place. 

Thus, the number of wounded constituents in a nucleon that passed through the 
nucleus of mass number A at a given impact parameter b is 

V--* „,,.> 2-|l-p,G(6)|-^-|l-p,G(6)l-< 

5--=- "'=^'<''> = 1 - [1 - Gibt ■ * ' 

i.e. the second term in Eq. (8). 
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